Model study on the photoassociation of a pair of trapped atoms 
into an ultralong-range molecule 
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Using the method of quantum-defect theory, we calculate the ultralong-range molecular vibra- 
tional states near the dissociation threshold of a diatomic molecular potential which asymptotically 
varies as —l/R 3 . The properties of these states are of considerable interest as they can be formed 
by photoassociation (PA) of two ground state atoms. The Franck-Condon overlap integrals between 
the harmonically trapped atom-pair states and the ultralong-range molecular vibrational states are 
estimated and compared with their values for a pair of untrapped free atoms in the low-energy 
scattering state. We find that the binding between a pair of ground-state atoms by a harmonic trap 
has significant effect on the Franck-Condon integrals and thus can be used to influence PA. Trap- 
induced binding between two ground-state atoms may facilitate coherent PA dynamics between the 
two atoms and the photoassociated diatomic molecule. 

PACS numbers: PACS numbers: 32.80.Pj, 33.80.-b, 34.10.-fx, 42.50.Ct 



INTRODUCTION 

Photoassociation (PA) spectroscopy Q] has developed 
into an important tool for studying properties of cold 
atoms and diatomic molecules. The formation of ex- 
cited diatomic molecular states through cold atom PA 
had been reported a long time ago Over the last 
decade, theoretical techniques for calculating PA rate co- 
efficients have been developed by many groups js] . Pho- 
toassociation spectroscopy enables a precise calibration 
of long-range interatomic forces. It leads to the forma- 
tion of translationally cold molecules and thus provides 
one possible route of creating a molecular Bose-Einstein 
condensate (BEC). Two-photon resonant Raman PA cou- 
pling can also be utilized for generating atom-molecule 
coherence 0, Q and massive entangled atoms 0. In 
fact, many-body atom-molecule coherence was recently 
observed |j| through the tuning of a Feshbach resonance 
(FR) @ caused by an applied magnetic field in a 85 Rb 
condensate. Quasibound molecular states can be formed 
by FR in ultracold atom-atom collision. FR has striking 
analogy with PA. In the first approximation, both are 
coupled two-state systems and mathematically identical. 
It is the magnetic field which couples the two states in the 
former case, while in the latter case it is the optical field. 
Both can be used to change the scattering length between 
two atoms 0. Two-color stimulated Raman adiabatic 
passage (STIRAP) Hill], [ll El via coherent PA using 
an appropriately tailored pair of laser pulses has been 
proposed for creating a molecular condensate from an 
atomic BEC. Stimulated PA of condensed atoms leads to 
"superchemistry" 0,E| where quantum statistics plays 
a crucial role in the production of molecules, and a super- 
position of macroscopic number of atoms and molecules 
|l6| . Although there has been no explicit experimental 
realization of a molecular BEC, coherent production of 



molecules in a single rotational-vibrational state has re- 
cently been experimentally demonstrated [t7| . The quan- 
tum dynamics (lif of coupled atomic and molecular con- 
densates reveals a lot of interesting physical effects, e.g., 
the recently observed frequency shifts contains the 
detailed information about the density of states of cold 
atoms in the quasicontinuum regime. An alternative ap- 
proach |2(J for creating a molecular BEC from atoms in 
a Mott insulating state inside an optical lattice was also 
proposed recently jU • 

In a PA event, two ground-state atoms, by means of 
a single-photon excitation during the collision, combine 
to form an excited diatomic molecule which can decay 
back to two atoms or to a different ground molecular 
state. Near-zero-energy vibrational states close to the 
dissociation threshold of an excited diatomic molecular 
level can be formed by PA of cold atoms. These states 
can extend to several hundred nanometers. PA absorp- 
tion line strength and linewidth strongly depend on the 
overlap integral, known as Franck-Condon (FC) integral, 
of these states with two atoms in the ground state. For 
these ultralong-range states, the "reflection approxima- 
tion" which is often used in molecular spectrum calcu- 
lations, breaks down. The reflection approximation re- 
places the actual integral by an integral over a 5 func- 
tion. This drastically simplifies the calculation by ruling 
out the necessity of calculating the wave function in the 
entire range. Instead, it only requires to calculate the 
wave function at a particular interatomic separation (the 
turning, or Condon point). This indeed seems to be a 
valid approximation for deeply bound molecular states 
with energy several tens of GHz below the dissociation 
threshold. For bound states with energies of the order of 
MHz or even kHz below the dissociation threshold, how- 
ever, this approximation does not seem to work. Calcu- 
lation of the actual Franck-Condon overlap integral is, 
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therefore, inevitable for such states. This presents a sig- 
nificant challenge as the standard numerical techniques 
for calculating the wave functions of these states are not 
known to be reliable. 

In this paper, we discuss a scheme for estimating the 
qualitative dependence of the required FC integrals. By 
using the quantum-defect theory (QDT) to find the ex- 
act solutions of molecular states in a long-range excited 
potential (-1/i? 3 ), we obtain near-zero energy (< MHz) 
vibrational wave functions that can extend from several 
hundred nanometers to micrometers. We employ these 
wave functions to evaluate the FC overlap integral in two 
situations: (1) for two atoms trapped in a harmonic po- 
tential and (2) for two free atoms. We analyze the effect 
of the trapping potential and the atom scattering length 
on the FC integral. We also evaluate spontaneous emis- 
sion linewidths of the excited molecular states. 

This paper is organized as follows. In the following 
section, we describe the model system and the basis for 
our study. Then in Sec. Ill, we discuss the relative mo- 
tional state of two ground-state (S+S) atoms trapped 
in a harmonic potential. In Sec. IV, we describe the 
method of QDT for calculating the vibrational states 
that can be formed by PA in the excited electronic po- 
tential which asymptotically corresponds to two (S+P) 
separated atoms. Photoassociation of two atoms inside 
an isotropic harmonic trap, with an emphasis on the use- 
fulness of the trap- induced two-atom bound states in one- 
and two-color coherent PA, is discussed in Sec. V. Wc 
present and discuss the results of our study in Sec. VI. 
We conclude in Sec. VII. 



THE MODEL SYSTEM 

Using the Born-Oppenhcimcr approximation, the 
Hamiltonian of our model with two atoms A and B can 



be expressed as H = H, 
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V 2 iel + V(R). (1) 



M and \x are the total and the reduced mass of the two 
atoms, respectively. The Laplacians V^ M and V^ ei , re- 
spectively, correspond to the center of mass and the rel- 
ative nuclear coordinate, of the two atoms. V(R) is the 
atom-atom interaction potential which is approximated 
as isotropic, i.e., as a function of the relative nuclear co- 
ordinate R = |Ra — Rb| only. H e \ is the electronic part 
of the total Hamiltonian and Vtrap = J2 a =A b ^trap with 
V t " ap as the trapping potential for atom a. The atom- 
field interaction is given by the standard dipole approxi- 
mation 
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where E represents the electric field of the applied laser 
field and d Q denotes the electronic dipole moment of 
atom a. In the limit of a low-intensity laser, this atom- 
field interaction can be considered as a perturbation. 

We choose molecular states as our basis functions. Ini- 
tially, the two atoms A and B arc in their ground s elec- 
tronic states. Thus their initial molecular state nota- 
tion is = \ns,ns, 2S+1 A g , u ,ei, J h Mi)\$i{R)), where 
A is the projection of the total electronic orbital angular 
momentum on the molecular axis and e$ is the relative 
energy of the two atoms. \^i(R)) is the initial relative 
nuclear wave function. The final molecular state can be 
represented by |$/) = \ns, np, Q, u>g , e v , J/, M f )\^ v (R)), 
where v denotes the vibrational quantum number, the 
energy E v of the vibrational state v lies at an energy 
e v = D f — E v below the dissociation threshold Dj. Here 
£l u ,g refers to Hund's case c of the interatomic potential; 
Ji(f) is the total angular momentum (consisting of elec- 
tronic orbital, spin, and molecular orbital angular mo- 
menta) of the initial (final) state with respect to a labora- 
tory axis and Muf\ is its projection on that axis. \^> V {R)) 
is the relative vibrational nuclear wave function. These 
basis states are similar to those used in Ref. [ij^. In 
this representation of the basis states, we have neglected 
the hyperfine interaction for the sake of simplicity. The 
dipole interaction term (J2J in this molecular basis, in the 
center-of-mass frame of the two atoms, then takes the 
form 
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where D(R) is the molecular dipole matrix element be- 
tween the molecular states A 9iU and £l u ,g, k^, and e are 
the wave vector and the polarization unit vector, respec- 
tively, of the laser field; and fx is the unit vector of the 
dipole moment. For two homonuclear alkali metal atoms, 
D(R) has been calculated by Marinescu and Dalgarno 
[23j. To leading order, it asymptotically varies as 1/R? in 
the form of D(R) = D Q + Di/R a . For the 351/2 - 3P 3/2 
transition in Na 2 , in atomic units Dq = —3.5007 and 
Di = 142.13 [13. Therefore, for R larger than 100a in 
unit of Bohr radius eto, the R dependence of D (R) of Na 2 
can safely be neglected and D(R) can be approximated 
as a constant Dq. All ultralong-range excited vibrational 
states of Na 2 we report here have negligible amplitude in 
R < 100; instead, their major amplitude lies at a sepa- 
ration much greater than a 1000 Bohr radii. 

To illustrate our study with concrete results, we specif- 
ically consider the transition 3 S U — > 0~ of Na 2 . The 
initial and the final internal states are assumed to be 
\ns, ns, 3 S u , ej, 0, 0) and \ns, np, 0~, e v , 1, 0), respectively. 
The potential of the 0~ and l u states of Na 2 has been 
explicitly tabulated by Stwalley et al. |24|. These states 
have their equilibrium positions at a separation larger 
than 60eto. Compared to the size of a ground diatomic 
molecule whose equilibrium position ranges in general be- 
tween 2 and 10 Bohr radius, these 0~ and X u molecular 
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states are truly long ranged. The PA spectrum for some 
vibrational states of ~ o f Na2 have already been ob- 
served experimentally |25j . 



TWO-ATOM STATE IN A TRAP 

We can choose the initial state of two ground-state (S) 
atoms in two different ways depending on whether the 
two atoms are trapped or are moving freely. The corre- 
sponding PA spectrum for these two cases can become 
markedly different. For two free ground-state atoms col- 
liding under the interatomic interaction, their asymptotic 
relative motion scattering wave function is given by 



* f = \ I — sm[k(R-a sc )] 



(4) 



which is an energy-normalized scattering state, e = 
fi 2 fc 2 /(2/x) is the asymptotic collision energy. At low en- 
ergies, for \R — a sc \ <C k^ 1 this wave function can be 
approximated as a straight line. Care should be taken 
to calculate the free-bound FC factor of ultralong-range 
molecular states, since the range of the vibrational state 
can exceed even at very low energies. 

For two atoms in an isotropic harmonic trap, their 
wave function is separable into the center of mass and 
the relative motion |2fil |. Approximating the atom- 
atom interaction V(R) by a regularized contact poten- 
tial 47rS 2 a sc (5 3 (R)/(2/x) (with a sc as the s-wave scatter- 
ing length) , the bound state of two atoms in a harmonic 
trap (of frequency ujt) was first derived analytically by 
Busch et al. |2^|. Let us define a characteristic length 
scale of an isotropic harmonic trap as at — y/h/JJwjt) 
and introduce a dimensionless quantity R — R/a t . Let 
f „, = RQ, H represent the relative motional bound-state 
wave function in the harmonic trap, with n t (which we 
will define below) being an integer quantum number re- 
lated to the three-dimensional isotropic harmonic trap. 
We can then express the bound state as 

$ nt = ^-V 2 Aexp(-R 2 /2)T(-v)U(-v, I R 2 ), (5) 

where A is a normalization constant having the dimen- 
sion of inverse of square root of volume and 



TABLE I: The energy eigenvalues E v of the lowest three har- 
monic trap-bound s-wave vibrational atom-pair states and 
their outer turning point R t for two values of harmonic fre- 
quencies. 

uj t = 2tt x 10 kHz ut = 2tt x 100 kHz 
n t {v) E v (MHz) R t (nm) E v (MHz) R t (nm) 



1 

2 



x 10 



T 



0.095 


515.44 


0.96 


164.25 


0.221 


786.35 


2.23 


249.89 


0.347 


985.28 


3.49 


312.78 
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FIG. 1: Two 23 Na atom s-wave bound-state wave function 
(in units of nm" 1 ' 2 ) in a harmonic trap with uj t = 2n x 100 
(kHz) (a) and uit = 2ir x 10 (kHz) (b). The individual curves 
are specified by the trap quantum number nt- 



where x = e/(hu} t )- We here define a dimensionless quan- 
tity £ s = a sc /a t . If £ s <C 1, then the eigenenergies can be 
approximated as 



+ 2n t + y/2/^^C^ +1/2 



hu t , n t = 0,1, 2,. ...(8) 



2hu t 4' 



(6) 



an effective quantum number for the relative motional 
eigenstate. [/(., ., .) is the confluent hypergeometric func- 
tion. The energy spectrum for such trap-bound two-atom 
states has been analyzed in detail before [2(| . The energy 
eigenvalue is given by the roots of the equation 



1 



F(-x/2 + 3/4) 

r(-a;/2 + l/4) " 72^ 



(7) 



Selected numerical results are presented in Tableland 
Fig. ^ Recently, the energy spectrum of two atoms 
in a harmonic trap has also been analyzed in detail in 
two papers |27| where it was shown that the validity of 
the solutions found in Ref. 26] critically depends on a 
characteristic length scale of the interatomic potential. 
This length scale is defined by Pe — (2fiCe/h 2 ) 1 / 4 , where 
C*6 is van der Waal's coefficient and /i = m/2 is the 
reduced mass of the two atoms. If both (3q at and 
a sc a t , then the solutions of Ref. |26j are physically 
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valid. If these validity conditions are not satisfied, then 
the use of an energy-dependent scattering length in the 
pseudopotential improves the solution as demonstrated 
in Ref . [27| . For a trap frequency uj t = 27r x 100 kHz 
and for 23 Na atoms, we find £ s = a sc /a t — 0.042 and 
fis/at = 0.13. Thus for our model system pursued in this 
study as described in Sec. II., we can safely use the so- 
lutions of Ref. 26] for very strong traps with frequency 
Lj t < 2tt x 100 (kHz). 



BOUND STATES IN THE -C 3 /R 3 POTENTIAL 

In order to compute the PA spectrum, one first needs 
to find the vibrational wave functions of the excited di- 
atomic molecule. It turns out, however, that it is not 
always necessary to find vibrational states for the entire 
range-the knowledge of wave functions near the vicinity 
of an outer turning point sometimes allows one to make 
a good estimate of the FC factor. In such cases, the FC 
factor can be approximated as the product of the two 
vibrational amplitudes at the outer turning point. This 
approximation, known as reflection approximation, relies 
on the fact that the molecule has its major probability of 
existence at the semiclassical outer turning point. Bar- 
ring the ground state, this is true for most of the low-lying 
vibrational states of a molecule. But, for the vibrational 
states near dissociation threshold, the probability can be 
distributed over a wide range. Therefore, for those states, 
the validity of this approximation becomes questionable. 
This was examined in the past by many authors for the 
case of bound- free transitions [2^] that defy the reflection 
approximation. 

In this study, we concentrate on near-threshold vibra- 
tional bound states in a potential that varies asymptot- 
ically as — 1/R 3 , which corresponds to a S+P excited 
molecular state between two atoms in the separated atom 
limit. Recently, these states have become very important 
in view of their wide accessibility in various PA processes. 
It is widely appreciated that there exists a class of pure 
long-range excited molecules whose equilibrium position 
is at a large separation [24|. The binding energies of 
these vibrational states can range from a few (Hz) to sev- 
eral hundred GHz. The lower the binding energy is, the 
longer the spatial extension of the states is. In calculat- 
ing the FC factor for these states with respect to the state 
of two ground state atoms moving either freely or in a 
trap, use of the reflection approximation becomes ques- 
tionable, particularly so for the vibrational states with 
sub- GHz energies. Therefore, it is necessary to calculate 
the vibrational states in the entire range or at least in a 
broader range near the outer turning point. 

In practice, it is difficult to accurately compute wave 
functions for vibrational states just below the dissocia- 
tion threshold by the standard numerical integration pro- 
cedure such as the Numerov method. Therefore, alter- 



native procedures such as the semi-analytic approach of- 
fered by the method of QDT need to be explored. While 
QDT has been successfully employed over the years for 
computing Rydberg-atom-like energy spectrum, its ap- 
plication to the problem of PA has not yet been fully 
explored. For Rydberg-atom-like systems, the relevant 
potential is the Coulomb potential that admits analytic 
solutions known as Coulomb functions which forms the 
basis for the application of QDT. For the excited state 
in the single-photon PA process as discussed here, the 
potential is of an attractive —1/R 3 type. Fortunately, a 
semianalytic QDT |29j of PA has recently been developed 
which is based on the semiclassical WKB-type solutions 
of a —1/R 3 potential. Exact quantum solutions of this 
potential have also been found recently |3fll3l| . For some 
diatomic molecules, the spacings of near zero-energy vi- 
brational states can be as small as a few kHz or even a 
few Hz, which leads to the interesting possibility that the 
free-bound PA linewidth of these states can exceed the 
spacings by several times even at a temperature as low 
as 1 /iK. Such a situation prevents the resolution of indi- 
vidual spectral lines. One naive suggestion to overcome 
this problem is to place two atoms inside a single trap, 
and to photoassociate them for different trap frequencies. 
Although PA in a trap may lead to a diffused continuous 
absorption spectrum, as we will show below, the use of 
different trap frequencies results in shifts of the peak of 
the spectrum. The energies of the different vibrational 
levels can be estimated from these shifts. 

We now briefly outline the use of QDT for computing 
near-threshold ultralong-range vibrational states. The 
exact scattering solution of a repulsive 1/R 3 was first 
obtained by Gao [22j , and the exact analytic solutions of 
an attractive 1/R 3 have been applied to a QDT formu- 
lation of diatomic molecular excited vibrational energy 
calculations before 30] . The scattering solutions of an 
attractive 1/R 3 potential has also been used before for 
a multichannel QDT formulation of atom-atom interac- 
tions in the presence of a dc electric field [13] • 

We next concentrate on the bound solutions of a 
— 1/R 3 potential. The asymptotic form of linearly 
independent base functions /; and gi satisfying the 
Schrodinger equation with this potential for energy E < 
can be expressed as 
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where Ws are chosen to be real functions 
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FIG. 2: The last 7 (from 33 to 39) s-wave vibrational wave 
functions (in nm" 1 ' 2 ) of the OJ state of Na2. The numbers 
assigned to the curves refer to the vibrational quantum num- 
bers. The range of the last bound state (v — 39) shown in the 
inset is so long that it exceeds /im by two orders of magnitude. 

with 9 = irv/2 + 7r/4 and 

oo 

£=F = X! b m exp(±imir/2), (11) 

m— — oo 
oo 

Tfe= = Yj (-l) m b m exp{Timir/2). (12) 

m— — oo 

Following standard conventions, we take S- 2 K 2 /(2/i) = 
—E, the bound state energy. The other parameters 
C±{v) 1 D(v), G{v) and b m are as defined in appendix 
A of Ref.|3lj. The condition for bound states is given by 

X -K° = 0, (13) 

where K° is the short-range K matrix determined by 
matching the solutions of short- and long-range poten- 
tials at a judiciously chosen separation, and x is given 
by 

X =-W f -/W g -. (14) 

With the help of QDT, the short range K matrix is ex- 
trapolated from the positive energy scattering solutions 
near the dissociation threshold. The bound-state wave 
function at discrete energy Ei can then be expressed as 

%(R) = f l -K°g l . (15) 

Selected numerical results for the ultralong range 
molecular states are given in Table ITT1 and Fig[3 

PHOTOASSOCIATION IN A HARMONIC TRAP 

A key factor of the PA spectrum strength is the FC 
integral. Ever since its introduction more than 70 years 



TABLE II: The energy E v of the last seven molecular vibra- 
tional levels of the 0g state of a Na2 for rotational quantum 
number I = and their outer turning point R t and the sepa- 
ration -R m ax at which the probability ^(i?)! 2 peaks. 

v Eg (kHz) R t (nm) i? max (nm) 



33 


1460 


162.3 


139.3 


34 


550 


224.8 


190.2 


35 


170 


332.3 


275.4 


36 


39.5 


540.3 


431.5 


37 


5.7 


1000 


794.0 


38 


0.3017 


2700 


1900 


39 


0.0003 


28200 


14200 



ago |3Jj | , the FC principle remains a paradigm in molecu- 
lar spectroscopy, particularly for the spectra of diatomic 
molecules. The essence for this principle relies on the 
wave mechanical property of molecules. The underly- 
ing idea is that when a molecule undergoes electronic 
transition from one state to another, the transition oc- 
curs much faster than the time scale of the vibrational 
motion of its constituent atomic nuclei. Therefore, just 
before and after the transition, the relative nuclear posi- 
tion remains almost the same. The wave mechanics of the 
molecule comes into the picture when one calculates the 
transition probability between the two states. In view of 
the FC principle, the strength of the transition is deter- 
mined by an integral over the overlap of initial and final 
wave functions of the nuclear motion. The line strength 
and the linewidth of molecular spectrum are proportional 
to this overlap integral. The ground vibrational state of a 
molecule has the maximum probability at its equilibrium 
position, while the excited vibrational states have their 
maximum probability near outer turning points. Accord- 
ingly, the relative location of the equilibrium positions is 
also important in determining which vibrational transi- 
tion would have a given intensity. 

For two atoms initially in a harmonic trap, we can 
define a Rabi frequency 

Q nt -v = h(nt\Va\v)\ (16) 
h 

between the trap-induced bound state \nt) of the two 
atoms and the excited molecular vibrational state \v) due 
to the interaction of the two-atom system with the PA 
laser field. The Rabi frequency Eq. (|16|l is proportional 
to the rate of coherent population transfer between the 
two-atom state \n t ) and the molecular state \v). Clearly, 
it is proportional to overlap integral 

/ *„*„ t cos(k L • R)dR, (17) 

which is the FC factor between the two states. Since 
both states \v) and \nt) are discrete (bound) and unit 
normalized, the FC factor in this case is dimensionless. 
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Here we wish to stress that the dynamics of photoasso- 
ciation of two trapped atoms is strikingly different from 
that of two free atoms. Since two trapped atoms form dis- 
crete relative motional states, there exists a well-defined 
Rabi frequency and the associated coherent dynamics. 
In contrast, for two freely moving atoms, the defini- 
tion of free-bound Rabi frequency is not so straightfor- 
ward, since Rabi frequency is usually defined between 
two discrete states. Instead, for low laser power, applying 
Fermi's golden rule, one can define free-bound stimulated 
transition rate 

Kr=Y l«f|e>| 2 , (18) 

which is proportional to the square of the free-bound FC 
factor 

Vv-e = J *„*e cos(k L • B)dR, (19) 

having the dimension of 1/s/e, i.e., |77„_ e | 2 has dimension 
of per unit energy. T^ -6 is the free-bound stimulated 
decay rate and does not characterize any coherent PA. 
In order to describe free-bound coherent PA dynamics, 
a quasicontinuum model [12J has been proposed. In this 
model, quasi-continuum is assumed to consist of many 
discrete equally spaced levels and thus enables to define 
a photoassociative Rabi frequency. Then the continuum 
limit is taken by allowing the energy spacing between 
the levels to approach zero. With this model analysis, 
it was shown |12| that the free-bound Rabi frequency is 
proportional to the square root of the collision energy. 
Thus in the zero energy limit, the free-bound Rabi fre- 
quency vanishes and so there is no coherent photoassocia- 
tive dynamics. Alternatively, the quasicontinuum model 
can be supplemented by introducing two-atom scat- 
tering states normalized within a sphere of arbitrary vol- 
ume V. After all the calculations are done, the contin- 
uum limit can be taken by allowing the volume V to go 
to infinity. It is thus argued that for a nondegenerate 
gas, free-bound Rabi frequency is much smaller than the 
corresponding bound-bound Rabi frequency. To transfer 
atoms into molecules by STIRAP, one of the necessary 
and sufficient conditions is the swapping of the bound- 
bound Rabi frequency into free-bound Rabi frequency in 
a counterintuitive way, i.e., for an initial duration of the 
STIRAP pulse, bound-bound Rabi frequency should be 
much larger than the free-bound one, and later it should 
be just the opposite. Therefore, according to the argu- 
ments put forward in Refs. [lolll^ . the STIRAP condi- 
tion in a nondegenerate gas would not be satisfied and 
so there is no STIRAP in a nondegenerate gas; though, 
there has been a debate [34| on the possibility of STI- 
RAP in a nondegenerate gas. However, it is also argued 
[lfij | that, for Bose condensed atoms and in the thermo- 
dynamic limit, the collective free-bound Rabi frequency 
becomes proportional to the number density of atoms 



and so it regains a finite (nonzero) value. Therefore, for 
condensed atoms, STIRAP becomes possible because of 
the symmetric state bosonic stimulation [l4| . In passing, 
we note that the possible effect of a trap on nondegen- 
erate free-bound-bound STIRAP was earlier pointed out 
in Ref.[l2. 

We emphasize that with trap-induced two-atom bound 
states, unlike in the case two free atoms, two-photon 
Raman-type photoassociative coherent dynamics is pos- 
sible. The conditions for STIRAP can also be satisfied 
for a sufficiently dilute nondegenerate atomic gas in a 
harmonic trap due to the existence of trap-induced two- 
atom bound states. However, the mean-field energy shift 
and inelastic collisions such as the three-body interaction 
may limit the efficiency of the STIRAP. 

RESULTS AND DISCUSSIONS 

Having discussed in the preceding section the advan- 
tages we may possibly derive from trap-induced atom- 
pair binding in performing two-color photoassociative 
STIRAP in a nondegenerate atomic gas to produce cold 
molecules, we now turn our attention to the effects of this 
pair binding on incoherent one-color photoassociative 
processes such as spontaneous and stimulated linewidths. 
The spontaneous linewidth for a transition v — ► n t is 
given by 

„,v — >n-t , .3 In |2|~, |2 i abound — bound { 

7sp 3^3 °l I*?"-"*! ' ^O) 

where u> vnt — uja — (^>v + w„, ) , uja is the frequency of 
the bare atom transition, Hu) v — E v is the energy of the 
bound molecular level v, hiu nt — E nt is the energy of the 
trap-bound level n t and 

7 bound-bou„d K ^uj 3 vv J(v\D(R)\v g )\ 2 (21) 

is the rate of spontaneous emission for transition from 
the vibrational state v in the excited molecular potential 
to the bound levels in the ground molecular potential. 
Lu VVg = uj a — {^v + w v ) , with Vg as the vibrational quan- 
tum number of the bound state \v g ) in the electronic 
ground molecular potential. For a bound-free transition 

7 ^ £ = 3^/ ^l^o| 2 |^-,| 2 & + 7, bound " bound (22) 

where uj ve = uja — (w» + e/K). The integration on the 
right-hand side of the above equation is over the distri- 
bution of collision energy e. 

The contribution to the total spontaneous line width 
due to transitions from the ultralong-range bound states 
in the excited molecular potential to bound states in 
the ground molecular potential is negligible. Because 
the turning points of even the least-bound state (longest 
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range) in a typical ground molecular potential that 
asymptotically becomes a van der Waals interaction of 
the form —Cq/R 6 are of the order of 1 nm, as can be es- 
timated from the binding energies of such states for L12 
and Na2 as calculated in Refs. [351I36T ]. In contrast, the 
excited ultralong-range vibrational states, as calculated 
here by the QDT method in Sec. IV, have their turn- 
ing points at a much larger distance of the order of 100 
nm. We further note from the plots of the wave func- 
tions of these ultralong-range vibrational states (Fig. [2J 
that, below 10 nm, the amplitude of these wave functions 
is vanishingly small. Therefore, the FC overlap integral 
of these excited ultralong-range states with the bound 
states in the ground electronic manifold would be negli- 
gible in comparison to that with the two-atom scattering 
states as well as trap-induced bound atom-pair states. 
We henceforth calculate the spontaneous emission width 
for a transition from a particular ultralong-range vibra- 
tional level in the excited molecular potential only to the 
trap-induced bound atom-pair states and to the two free- 
atom states (scattering states) while neglecting all other 
molecular bound-bound transitions. In this context, a 
pertinent question that may arise is whether such tight 
traps would have any influence on the molecular bound 
states in the ground molecular electronic manifold. As we 
have already mentioned, the least-bound or the nearest- 
to-zero-energy molecular bound states supported by van 
der Waal's potential —Cq/R 6 have their outer turning 
point at a typical separation of 1 nm. If we now esti- 
mate, for a typical value of Cg = 1500 a.u., the van der 
Waal potential energy at R = 1 nm, and compare this 
with the relative harmonic trap energy of two atoms at 
the same separation in an isotropic harmonic trap with 
harmonic frequency 100 kHz, we find that this trapping 
potential energy is smaller than the van der Waal po- 
tential by at least two orders of magnitude. So, for har- 
monic trapping frequency lot < 2tt x 100 kHz, the trap 
has negligible influence on the molecular bound states in 
the ground electronic molecular potential. This is why 
in writing the expressions for 7^™* and yj£* e in Eqs. 
(|20|l and l|22|) . respectively, we have retained the same 

7 bound~bound m both express i ons . 

In Sec. Ill, we have described two-atom bound states 
ty nt (R) formed inside a harmonic trap. Unlike the free- 
bound case, the FC factor for two trapped atoms can 
be controlled by tuning the trap parameters. The lower 
the trap frequency is, the longer the extension of the 
trap-induced two-atom bound state. Thus, by lower- 
ing the trapping frequency, the FC factor for the higher 
molecular vibrational states can be increased. Therefore, 
the trap-induced bound state can facilitate probing of 
ultralong-range molecular states. Figure ^ exhibits the 
lowest three ^ nt for two values of trapping frequency 
uj t = 27r x 100 or 2n x 10 kHz. Since £ s -C 1, they de- 
viate very little from the harmonic-oscillator states. The 
bound-state energies and the outer turning point of these 
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FIG. 3: Dimensionless Franck-Condon integral rj vnt be- 
tween molecular bound v = 35 and the trap-bound rit = 
(ground) states as a function of positive scattering length a sc 
in nanometers for a trap frequency of uit = 2n x 10 kHz. The 
inset shows the variation of trap-bound ground-state energy 
eo (2-kvq) as a function of a sc . 

states are given in Table U 

The wave functions (R) of the last seven vibrational 
states from v = 33 to v = 39 calculated by the QDT 
method as outlined in Sec. IV are plotted in Fig. [21 
Their vibrational energies E v , outer turning points Rt, 
and the internuclear distances i? max at which \^f v \ at- 
tains maximum are given in Table [HJ To evaluate E v , 
we employ Eq. I|15fl with a constant K° = 0.623. The 
vibrational energies of these states have also been calcu- 
lated previously by a similar method |30|. From Table ITT1 
it is clear that the molecule in these states would spend 
most of its time at an internuclear distance larger than 
100 nm. In terms of typical molecular scale, this length 
scale is enormous. Now, a question naturally arises: is it 
possible to probe these states by PA? To seek an answer 
to this question, we carry out an explicit calculation of 
the FC factors in two different cases: (1) both atoms be- 
ing trapped in a single harmonic potential and (2) the 
two atoms are in a free scattering state. 

In Table III, we tabulate the values of the square of 
the FC factor for selected states. Why and how much 
these values differ from one pair of states to another can 
be naively understood from the positions of their outer 
turning points and the shape of the respective wave func- 
tions. For example, from Tables III and ITT1 we note that 
the states ^ nt =o [f° r — 2-zr x 100 kHz] and 'fu^s 
have their outer turning points at a comparable distance. 
Therefore, these two states should have maximum over- 
lap. 

In Figs. and we display the variations of the 
bound-bound FC factor as a function of the positive 
and negative scattering length, respectively. Although 
as discussed earlier that for a large positive a sc , energy- 
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FIG. 4: The same as in Fi.g[3]but for a negative scattering 
length. 



TABLE III: Square of Franck-Condon integral \r) vnt | 2 between 
molecular bound v and trap-bound states nt . 

Lut = 2n x 10 (kHz) u t = 2tt x 100 (kHz) 

v nt = nt = 1 nt = 2 nt = nt = 1 nt = 2 

33 0.075 0.074 0.059 0.221 0.122 0.009 

34 0.102 0.080 0.039 0.024 0.066 0.019 

35 0.116 0.004 0.033 0.000 0.214 0.021 

36 0.009 0.022 0.000 0.000 0.013 0.093 



dependent pseudopotential should be used in the calcu- 
lation of trap-induced atom-atom bound states, we sim- 
ply took the energy-independent pseudopotential for all 
the calculations. This approximation underestimates the 
actual eigenenergies and hence underestimates the ac- 
tual range of those states. Therefore, for large positive 
a sc , the actual value of the FC factor should be some- 
what larger. However, we took this approximation for 
our model calculation and comparison because the over- 
all pattern of variation of the FC factor as a function of 
a sc should remain the same. 

When two atoms are not trapped, the FC factor be- 
comes strikingly different. Since the FC factor in the 
two cases have different dimensionality, it is not easy to 
have a direct comparison between the two. However, we 
can make a comparison by comparing their spontaneous 
linewidth which is proportional to the respective FC fac- 
tor. For the sake of this comparison with the bound- 
bound v — ► n t transition, we have assumed that the the 
velocity of the free atoms follow a Maxwellian distribu- 
tion at a temperature T = E nt /kB in calculating the cor- 
responding free-bound transitions, where Ub is the Boltz- 
mann constant. The results are tabulated in Table llVl 
The free-bound FC factor of Eq. (|19|l depends critically 
on the locations of the antinodes of the sinusoidal scat- 
tering state of the two atoms relative to the outer turning 



TABLE IV: Spontaneous linewidths 7^T* n * for transitions 
from excited molecular bound states v to trap-bound states nt 
are compared with corresponding linewidth 7sp"* e for molecu- 
lar bound-free transition. The frequency of the harmonic trap 
is Lu t = 2tt x 100 kHz. 

"7(2tt) fc 7 s V7(27r) in (kHz), k B T = htu nt 



n t = 


= 


nt = 


= 1 


n t 


= 2 


v v — » nt 


v — > e 


v — > nt 


v — > e 


v — > n t 


v — > e 


33 4278 


28985 


2883 


40422 


2359 


47139 


34 464 


26591 


5532 


37216 


1282 


43467 


35 0.002 


23248 


276 


30781 


4150 


34831 


36 0.011 


6449 


40.0 


9434 


245 


11007 


37 0.001 


3429 


6.3 


4106 


29.6 


4249 


38 0.000 


599 


0.5 


623.7 


2.4 


735.2 


39 0.000 


10.2 


0.002 


7.8 


0.009 


7.4 



point of the molecular bound state. If the first antinode 
is located near the turning point, then the free-bound FC 
factor would be greater than if it is located far away from 
the turning point. 

From Table |IV| we note that the tabulated bound- 
bound spontaneous linewidths are much smaller than 
their respective free-bound values. This suggests that if 
we PA two trapped atoms, and if the resulting ultralong- 
range excited molecule remains trapped under the same 
or a different trapping potential, then the spectral line 
associated with the formation of the molecule may be 
resolved. 

The trap frequencies we have considered for our nu- 
merical illustration correspond to very tight traps. In 
this section, so far we have restricted our treatment to 
an elementary PA process of only two atoms inside an 
isotropic harmonic trap. Now, the question obviously 
arises how the tight traps would affect the PA of an en- 
semble of atoms in a trap, in particular of an atomic BEC 
for its potential use in forming a molecular BEC by collec- 
tive photoassociative STIRAP. In this regard, it becomes 
necessary to consider various effects due to mean-field 
interactions and inelastic three-body recombination pro- 
cesses |ll|, |l7| . Let us first discuss the possible effect of a 
tight trap on the efficiency of STIRAP. The presence of a 
trap introduces an additional length scale in the system, 
namely, the size of the harmonic oscillator ground state 
a t . For a typical magnetic trap in use today with Rb or 
Na atoms, trap frequency ui is of the order of 2ir x 100 Hz. 
The frequencies we have considered here are 10-100 kHz, 
usually associated with strong far off resonant optical 
dipole traps. A condensate in such traps would typically 
have dimensions smaller than an order of magnitude of 
that in a magnetic trap. One can then argue that, for a 
typical number of atoms N ~ 10 6 , the number density of 
atoms in such a tight trap would also be higher by two 
to three orders of magnitude. Therefore, the strength 
of the mean-field interactions between atoms and also 
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between atoms and STIRAP-produced molecules will be 
greatly enhanced. This increased mean-field interaction 
will lead to large frequency shifts and broadening in two- 
photon STIRAP PA spectrum, it also leads to a dramat- 
ically increased loss of atoms or molecules due to density 
dependent three-body collisions. Therefore, as already 
suggested in Ref. [37j , our study applies to cases of mod- 
est atom number densities associated with condensates of 
smaller numbers of atoms, of the order of 1000 or fewer. 

The enhancement of the Franck-Condon integral due 
to trap-induced binding between two ground-state atoms 
depends on the location of the outer turning point in the 
excited (S + P) molecular potential relative to the value 
of at- Our calculations reveal that the amplitude of the 
ground state of the two-atom (S-S) bound state attains a 
maximum at about at/2. Therefore, the Franck-Condon 
integral will be maximal for transitions to states in the 
excited potential whose outer turning points lie close to 
at/2. For a typical magnetic trap of 100 Hz, a t is about 
1 /jm, while the outer turnin g p oint of the excited vi- 
brational state used in Refs. 1 171 is at about 3 ran. 
The finding of Drummond et al. |ll| that a weak trap 
reduces the efficiency of STIRAP, we believe therefore, is 
due to this huge mismatch between a t and the outer turn- 
ing point. They have used the excited vibrational state 
whose free-bound transition frequency is 23 cm -1 below 
the dissociation threshold. This state is far-off the dis- 
sociation threshold in comparison to the ultralong-range 
states we have described in the paper. To achieve a signif- 
icant enhancement, the free-bound transition should in- 
volve the highest lying excited vibrational states, so that 
their outer turning points lie close to a t /2. Normally, 
the excited ultra-long range vibrational states have their 
outer turning points at about tens to several hundred 
nanometers. In order to have large FC overlaps between 
these states and the trap-induced two-atom bound states, 
the trap frequencies should be much larger than 100 Hz, 
i.e., we should use very tight traps. As long as the con- 
dition (ndg C <C 1) remains fulfilled, the mean-field effects 
and three-body recombinations do not lead to substan- 
tial degradation of the efficiency of the STIRAP, as the 
gas remains sufficiently dilute and the interactions would 
be mainly of a two-body type, and only a small frac- 
tion of atoms will be pair correlated with trap-induced 
binding, they will experience enhanced FC coupling in 
a STIRAP. Presently, the nature of collision interaction 
between photoassociated excited or STIRAP-produced 
ground molecules and individual atoms is not precisely 
known. Therefore, we will omit a detailed study of how 
three-body recombination will affect the efficiency of the 
STIRAP. Intuitively it can be understood that with a 
high number density, the increased three-body recombi- 
nation will lead to increased inelastic losses, which coun- 
teract any enhancement effect by tight trap. However, it 
can also be intuitively argued that for a very low density 
of atoms and molecules, such three-body recombination 



losses may perhaps be overcome in achieving a consider- 
ably efficient STIRAP. 

In view of the foregoing discussion, it now appears that 
the best way out to avoid mean-field shifts and three- 
body effects is to photoassociate atoms in a Mott insu- 
lating state of atoms in an optical lattice |2l|. In such 
a structured state of atoms, a pair of atoms can be ar- 
ranged to occupy a single site of an optical lattice. 

Before concluding this section, we would like to dis- 
cuss another important application of the ultralong-range 
excited vibrational states in coherent laser spectroscopy 
of yet a different type of potentially ultralong-range 
molecules, namely, molecules that can be formed in a 
controlled crossing of a Feshbach resonance 0- These 
exotic molecules are produced in magnetic-field-induced 
Feshbach resonance by the nonadiabatic coupling be- 
tween two close molecular potential curves, one of which 
acts as an energetically closed channel while the other as 
an open asymptotic collision channel. Coherent trans- 
fer of atoms into a molecular vibrational state close to 
the open-channel threshold occurs through the tuning 
of a magnetic field. Recently, these vibrational states 
have been numerically calculated [13 and shown to have 
an extent ranging from many tens to many hundreds 
of nanometers. Although, these molecules are transla- 
tionally cold, they are not generally cold in their vibra- 
tional states. Obviously, any two-photon resonant coher- 
ent Raman spectroscopy of these molecules would involve 
higher-lying excited vibrational states as we have calcu- 
lated by the QDT method. Thus, in addition to our find- 
ing of the trap enhancement in coherent PA, our study of 
the QDT approach in finding ultralong-range wave func- 
tions is important for applying the STIRAP to the newly 
attained capability of a controlled crossing of a Feshbach 
resonance, and the attempt to produce translationally as 
well as vibrationally cold molecules. 

CONCLUSIONS 

We have developed a scheme for estimating the FC 
factor involved in a single-photon PA of two atoms into 
an ultralong-range molecular state. In this approach, 
we have made two simplifications: first, the initial state 
of the atoms are based only on the asymptotic form at 
a large interatomic separation, which in the low-energy 
limit, only depends on the s-wave scattering length and 
second, for the excited molecular state (corresponding 
to a S+P asymptote), we have used the exact solutions 
of the asymptotic long-range molecular potential —1/R 3 , 
taking advantage of the QDT formulation. We have cal- 
culated the wave functions of the molecular states that 
are very close to dissociation threshold. We have cal- 
culated the FC overlap integrals between these molec- 
ular states and the initial states of the two ground- 
state atoms. We have shown that when two atoms are 
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trapped in a single harmonic potential well, the trap fre- 
quency can significantly affect the integral and hence the 
PA rate. We suggest that the ultralong-range excited 
molecule can be formed by PA of two ground-state atoms 
inside a trap and simultaneously switching on another 
trap that can confine photoassociated molecules. In this 
way, the damping of the excited molecule can be min- 
imized. We have also demonstrated the dependence of 
the FC integral on both the sign and the magnitude of 
the atomic scattering length. To form a ground-state 
molecule by stimulated Raman PA, the excited vibra- 
tional state should have sufficient amplitude in the short- 
range region such that it can possess a significant over- 
lap with the ground molecular vibrational state. For the 
model system studied in this paper, this seems practi- 
cally impossible. However, since trap potentials have a 
significant effect on the overlap integral as well, it can 
perhaps also be adjusted to affect stimulated Raman PA. 
In particular, when performing coherent PA in a nondc- 
generate gas, we have shown that a tight trap does of- 
fer an advantageous situation in comparison to atoms in 
free-scattering state. However, as discussed in the pre- 
ceding section, the mean-field shifts and three body re- 
combination effects may ultimately limit the efficiency 
for employing the STIRAP process to form a molecu- 
lar BEC, and it appears that the best way to exploit 
the trap-induced enhancement of coherent coupling for 
producing a molecular BEC by STIRAP with no ad- 
verse effects from three-body recombination and mean 
field shifts is to photoassociate atoms in the Mott in- 
sulating state with more than one atom per well of an 
optical lattice as first investigated in Ref. |20||. To this 
end, we would also need an accurate method for calcu- 
lating the ultralong-range excited molecular vibrational 
states. Very recently, far-off-resonance trap for atoms 
inside a cavity in the strong-coupling regime has been 
experimentally demonstrated j^a]. This trap is insensi- 
tive to the internal atomic states and has strong or tight 
trapping capability. Only a few number of atoms can 
be captured in such a trap for a considerable duration. 
Trap-induced enhancement may be experimentally stud- 
ied in photoassociation inside such a tight trap without 
any side effects of mean-field interaction and three-body 
recombination. Apart from the usual PA spectroscopy, 
our study can also find applications in the future studies 
of coherent laser spectroscopy of the recently discovered 
ultralong-range Feshbach molecules Q. 
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